Nucleosynthetic Layers in the Shocked Ejecta of Cassiopeia A 



Karl Isensee 1 , Greg Olmschenk 1 , Lawrence Rudnick 1 , Tracey DeLaney 2 , Jeonghee Rho 3 , 
J.D. Smith 4 , William T. Reach 5 , Takashi Kozasa 6 , and Haley Gomez 7 

ABSTRACT 

We present a 3-dimensional analysis of the supernova remnant Cassiopeia A 
using high resolution spectra from the Spitzer Space Telescope. We observe su- 
pernova ejecta both immediately before and during the shock-ejecta interaction. 
We determine that the reverse shock of the remnant is spherical to within 7%, 
although the center of this sphere is offset from the geometric center of the rem- 
nant by 810 km s _1 . We determine that the velocity width of the nucleosynthetic 
layers is ~1000 km s _1 over 4000 square arcsecond regions, although the velocity 
width of a layer along any individual line of sight is <250 km s _1 . Si and O, which 
come from different nucleosynthetic layers in the progenitor star, are observed to 
be coincident in velocity space in some directions, but segregated by up to ~500 
km s _1 in other directions. We compare these observations of the nucleosynthetic 
layers to predictions from supernova explosion models in an attempt to constrain 
such models. Finally, we observe small-scale, corrugated velocity structures that 
are likely caused during the supernova explosion itself, rather than hundreds of 
years later by dynamical instabilities at the remnant's reverse shock. 
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1. Introduction 



The supernova remnant Cassiopeia A (Cas A) is a unique astrophy sical laboratory due 
to its young age (^340 years - iThorstensen et al.ll200ll ; iFesen et al.ll2006l ) and small distance 



(only 3.4 kpc - iReed et al.l Il995l ). The remnant is just entering it s Sedov- Taylor phase 



so emission from both forward and reverse shocks can be detected (IHughes et al. 1 120001 ) . 
Emission at most wavelengths, including most of the infrared, is dominated by a ~120" 
radius "Bright Ring". The Bright Ring is formed when supernova ejecta encounter Cas A's 



ejecta rich in 0, Si, S, Ne, Ar, Ca, and Fe 


Chevalier & Kirshner 


1978t Douvian et al. 


1999 


Hughes et al. 12000: Willineale et a 


.2003 


Hwane & Lamine 112003c Lamine & Hwane 


2003 


Morse et al. 2004: Ennis et al. 2006 


)• 



Studies of optical light echoes from the the explosion near peak light have led to the 
observation of weak hydrogen lines, indicating a supernova Type lib origin for Cas A 
(IKrause et al. 1 120081 ) . Cas A's progenitor was therefore a red supergiant that had lost 



most, b ut not all, of its hydr ogen envelope. X-ray studies indicate a total ejecta mass of 
2M^ ( Willin gale et al.l 120031) . If one adds to this the mass of the central compact object 



(jChakrabarty et al. 



20011 ) . Cas A's progenitor had a total mass of at least 4M immediately 



of 



(Youne et al. 


2006; 


Vink et al. 


1996) 



Although Cas A's appearance is dominated by recently shocked ejecta, it also con- 
tains emission that is not the result of collisional ionization at t he reverse shock, but pho 
toionization by UV and X-ray emission fro m the shocked ejecta ([Hamilton fe Sarazinlll984 
Hamilton fe Fesenlll998t ISmith et al.l 12009). This material i s seen toward the cen tral region 



of the remnant at low radio frequencies ( 



20081 : ISmith et al.ll200a IPeLanev et al.ll2010 



Isensee et al. 



Kassim et al.lll995) and in the infrared (IRho et al. 



20101 ) , and was demonstrated to 



be at lower densities and ionization state than recently shocked material on the Bright Rin g 
through a combination of Doppler analysis and line ratio measurements (ISmith et al.ll2009l ). 
These ejecta are often referred to as "unshocked ejecta" since they have yet to encounter the 
remnant's reverse shock. That is not an accurate label, since Cas A's forward shock and a 
reverse shock interacted with the ejecta during the supernova explosion itself. 
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1.1. Previous 3D Maps 



3D Maps of Cas A have been made in the optical, infrared, and X-ray. Doppler re- 



const ructions in the optical used S and emission lines ([Lawrence et al.l Il995l ; iReed et al. 



19951 ) and showed that ejecta on the Bright Ring lie on a roughly spherical shell but are 
not uniformly distributed on that shell - most of the ejecta lie nearly in the plane of the 
sky. They also observe that the center of the sphere is offset from the geometrical center 
of the spherical shell by ~0.36pc along our line of sight. This indicates that the ejecta are 
not traveling at the same velocity in all directions, which is c onsistent with previous results 



which indicated an asymmetric expansion for the ejecta (e.g. iBraun 1119871 ; IWillingale et al. 



20021 ). These 3D reconstructions give us a selective snapshot of ejecta because only material 
that has recently encountered the remnant's reverse shock will emit strongly in the optical. 



DeLaney et al.l (120101 ) create d a 3D infrared and X-ray map of Cas A from a Spitzer 



Space Telescope spectral cube3- Isensee et al. ( 2010 ) used a similar IR data set, but at 
higher spectral resolution, to make a 3D map of ejecta in the center of the remnant. The 
advantage of these IR maps lies in the fact that much of the ejecta in the IR will be detectable 
both before and after they interact with the reverse shock. Both studies found a similar 
distribution of ejecta to that seen in the optical where the center of expansion is offset from 
the geometrical center of the remnant both in projection and along the line of sight. These 
works were able to study the relationship of several nucleosynthetic layers and are discussed 
in the next section. 



1.2. Separation of Nucleosynthetic Layers 



Si and O emission are observed to be co-located in most regions (e.g. lEnnis et al.l 120061 ) 
in both the X-ray and infrared. This indicates that the two layers have very similar velocities 
(less than 80 km s" 1 differen ce). However, eviden ce of layer diff erentiation is foun d in some 
directio ns in the X-ray (e.g 



IR (e.g. IPeLanev et al.l 12010 



lushes et al. 



Isensee et al. 



2000), the optical dFesen et al.ll2006h . and the 



20101 ). which was likely caused by the different 



layers of the star being ejected at different velocities in those directions, thus encountering 
the remnant's reverse shock at different times. 

It should be emphasized that we can only observe mixing or separation in velocity space. 
We can easily detect any velocity gradients in the supernova explosion since we can detect 
Doppler velocities of <100 km s -1 in the IR, while typical observed velocities and velocities 



Movies showing this 3D structure are available at 



http: / / chandra.harvard.edu/ photo/2009 / casa2 / animations.html 
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predicted by models are an order of magnitude larger (e.g. lHammer et al. II2010I ). However, 
we cannot detect any initial spatial separation of the nucleosynthetic layers - simulations 
pred ict that the relevant n ucleosynthetic layers will be < 10 11 cm thick prior to the explosion 
(e.g. IJoggerst et al. 1 120091 ) . but the typical ejecta clump size of <1" corresponds to ~ 10 16 
cm at Cas A's distance. Therefore, we cannot differentiate between a situation where two 
nucleosynthetic layers were separated during the supernova explosion but ejected at the 
same velocity, and one where the two layers were completely mixed during the explosion 
and ejected at the same velocity. But, if we observe two layers that are currently separated 
in velocity space, we know that they were separated during the supernova explosion itself 
because, to the best of our knowledge, there is no mechanism that will impart different 
velocities to spatially overlapping elements. 



1.3. Geometrical Asymmetries 



Supernova explosion model s predict substantial asymmetries due to effects such as rota- 



tion as well as instabilities (e.g. iBlondin et al. 1120031 ; iBurrows et al. Il2007t lHammer et al. 



20101 ). Observations of both supernovae and supernova remnants have confirmed this picture. 
Spectropolarimetric observations of unresolved supernovae have shown that all observed core 
collapse supernovae c ontain intrinsic pola rization, indicating that there is a departure from 
spherical symmetry (IWheeler et al.l 120051 ) . Although an axis-symmetric geometry, proba- 
bly induced by jets, can be used to explain some features in some core collapse super- 
novae, significant depar tures from axial symmetry are needed to explain most observations 
rtWang fc Wheeler! l2008h . 



IR 3D maps of supernov a ejecta in Cas A have found major asymmetries, both on global 
scales ( DeLanev et al.l 2010 ) and for smaller subsets of ejecta in the the entire supernova 
remnant ( llsensee et al. 1120101 ). These asymmetries are not immediately apparent in the 
visual appearance alone because the highly spherical reverse shock creates a large selection 
effect in that we can only observe ejecta near the shock at most wavelengths. 

In this paper, we present an analysis of high spectral resolution Spitzer mappings of 
the ejecta on the B right Ring of Cas A. This data set is an extension of that used by 
Isensee et al. I ( 120101 ). and it contains regions with both recently shocked and interior ejecta. 
In §2 we present the observations and discuss the methods used in our analysis. We describe 
those results in §3 and discuss the physical implications in §4. 
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2. Observations and Analysis 

The Spitzer Infrared Spectrograph (IRS) was used on August 30, 2007 to create spectral 
maps of select relatively bright regions in Cas A whose locations are shown and labeled in 
Figure [TJ High-resolution spectra (R~600 for all wavelengths) were taken from 20-35/mi 
using the Long High (LH) module in all regions and from 10-35/mi using both the Long 
High and Short High (SH) modules in the Southwest region. The FWHM of unresolved 
spectral features in these observations is about 0.06/nn at 35/nn and about 0.02/im at 13/im. 
The LH data were taken using a 61 second exposure at each position while the SH data 
were taken using a 31 second exposure at each position. The pixel scale of the observations 
is ~1.25" and ~2.5" for the SH and LH modules respectively. The background, which was 
taken from 3 separate 61 second observations adjacent to the remnant, was subtracted and 
3D cubes were created usin g the S19 version of the IRS pipeline and the CUBISM software 



package (ISmith et al.l 120071 ). The uncertainties in flux for each line of sight were calculated 



from the IRS pipeline using standard error propagation of the BCD level errors. 

The undersampling of the IRS modules limits our uncertainties. This is a systematic 
error that exists in the wavelength calibration data themselves, and is worst at the short- 
wavelength end of both the SH and LH modules. Our obtainable wavelength accuracy is 
limited to roughly 1/2 of a spectral bin, or about 100 km s -1 . The relative wavelengths for a 
given line can be measured with higher accuracy from position to position or within multiple 
Doppler components in a given position. 



2.1. Spectra 

Cas A's infrared spectrum is dominated by bright emission lines as shown in Figure [2j 
The LH observations contain lines from [O IV] at 25.9yum, [S III] at 33.48/im, and [Si II] at 
34.81/xm. We tentatively identify the line near 23/mi as the 22.9/xm [Fe III] line. The lines 
observed in the LH module typically have peak fluxes from 100-10,000 MJy sr^ 1 , with an 
rms noise of ~20 MJy sr _1 . The SH observation contains lines from Ne, S, and Fe. Typical 
peak fluxes are ~300 MJy sr _1 . 



2.2. Doppler Deconvolution 



We performed a Doppler deconvolution of the spec tral lines for each line of sight from 
each ion separately using a spectral CLEAN algorithm (jPing et al. Ill999[ ). We determined 
the uncertainties in Doppler velocity for each Doppler component by applying the spectral 
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CLEAN to synthetic line data with a realistic range of signal to noise ratios, and using line 
free data in order to model the noise. For both procedures, we use identical techniques to 
those described in llsensee et al. I ( 120101 ). 



Using synthetic data, the uncertainty in velocity for a single, isolated Doppler component 
was determined to be <25 km s _1 , however, we could not differentiate two components from 
one another along the same line of sight that were within 65 km s _1 of one another. Therefore, 
uncertainties in the absolute velocities are limited by the systematic errors in the calibration 
of ~100 km s" 1 rather than random uncertainties. 



As in previous studies (e.g. iReed et al.lll995l ; iDeLaney et al.ll2010t llsensee et al. 1120101 ) 



we assume that the ejecta have been freely expanding at a constant veloc ity in order to 



determ ine their spatial coordinate perpendicular to the plane of the sky. IDeLaney et al. 



( 120101 ) demonstrate that this is a good assumption for IR emission by showing that nearly all 
ejecta plotted on a Velocity vs. Radius plot fall on a semi-circle. We make a similar Velocity 
vs. Radius plot in Figure [3] from our data set an d find that ejec t a in t hese regions fall on a 
semi-circle that is consistent with that found by IDeLaney et al.l ( 120101 ). The assumption of 
constant velocity is still valid despite the fact that the ejecta were likely decelerated during 
the supernova explosion itself. The subsequent behavior after the explosion is essentially 
identical to free expansion at a reduced velocity because any deceleration happened near 
t=0, z=0 where z is the spatial coordinate perpendicular to the plane of the sky. 



3. Results 
3.1. 3D Maps 

We plot the Doppler components from both the 25.89/im [O IV] and the 34.81/mi [Si II] 
lines in Figures HJ EJ and for the Southwest, Northeast, and Southeast regions respectively. 
We converted our velocity coor dinates to spat i al coo rdinates using the conversion factor 



between the two determined by IDeLaney et al.l (120101 ) . This conversion factor is more ap- 



propriate than one calculated from our own data since it uses data from the entire remnant 
rather than a few select regions. The flux from each component is displayed by varying the 
transparency; the brightest 3D pixel (or "voxel") for a given ionic line is 80% opaque, while 
the opacity of all other voxels is linearly scaled downwards as a function of the intensity of 
the Doppler component. We have not plotted very weak ejecta with total fluxes less than 
15% that of the brightest velocity component. The other strong line, the 33.48/xm [S III] 
line, is from the same nucleosynthetic layer as Si and traces out nearly identical structures 
to the 34.81//m [Si II] line. Therefore, we do not show it here. The low density of the interior 
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ejecta (~100 cm -3 , ISmith et al.l 120091 ) implies that self-absorption within the ejecta will be 
minimal along the line of sight. We note that we are likely only observing the densest ejecta 
material, whether it is shocked or interior ejecta, since the emissivity should scale roughly 
as the density squared in both cases. 

The ejecta in the SW region form a distinct shell-like structure. The O and Si ejecta 
all lie along the same shell, although they fill different parts of the shell. Averaging over 
the entire SW region, there are an averge of only 1.3 Doppler components per line of sight, 
indicating that the thickness of the shell is <250 km s _1 along any given line of sight. The 
brightest ejecta in the NE region also forms part of a shell, but there are substantially more 
dim ejecta inside of the shell than in the SW. O and Si lie both on the shell and inside of the 
shell, although it appears that some of the O and Si is systematically separated in velocity 
space (see § 3.6). The SE region consists of an irregularly shaped region of both Si and O 
emission in addition to a region dominated by O emission in the western most part of the 
region. 



3.2. Iron 

We observe [Fe II] with the SH module at 17.9/im in the SW. We plot the Doppler 
components from this line with Si and O emission as shown in Figure [7J The SH data were 
binned 2 by 2 pixels to increase the signal to noise ratio. It is clear that the Fe lies on the 
Si+O shell described in the previous section. 

We confirm that the 25.9/im line is the 25.89/iin [O IV] and not the 25.98/iin [Fe II] line 
by comparing the Doppler structure of the 25.9/iin to that of the 34.81/iin [Si II] line for 
several lines of sight in the SE, SW, and NE. As an example, we display the results for one 
line of sight with strong 17.9/im Fe emission in Figure We show the Doppler structure 
for the 25.9/mi line under the assumption that it is all [O IV] and all [Fe II]. We obtain an 
excellent match under the assumption of [O IV], but a poor match under the assumption of 
[Fe II] even along this line of sight where we see relatively strong Fe at other wavelengths. 
We find no ev i dence for Fe at 25.9/im for 10 other lines of sight, confirming the results of 
Isensee et al. ( 201ol ). Therefore, we assume for the remainder of this paper that the 25.9/um 



line is entirely due to [O IV] emission. 
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3.3. Comparison to X-ray Emission 



We compared the loca tions of the IR ejecta to X-ray ejecta detected in the 2004 Chandra 



observations of Cas A (see lHwang et al. II2004J ). The spectral resolution of the X-ray images 



is not sufficient to accurately determine the Doppler velocity for most lines. Therefore, we 
show the 3D location of the X-ray ejecta as planes perpendicular to the plane of the sky in 
Figure [9j These planes represent the forward edge of the X-ray ejecta as seen in Figure [TJ 
We know that the X-ray material has been recently shocked since the ejecta will only be at 
the appropriate ionization states and temperatures if it has recently encountered the reverse 
shock - the ejecta will ionize up to states that are difficult to observe in the X-ray within 
~100 years (TMazzotta et al. Ill998h . 



Most of the bright IR ejecta are immediately interior to the leading edge of the X-ray 
material. Since it takes some time for a plasma to up ionize to ions visible in the X-ray, this 
is consistent with the picture that the brightest IR ejecta have been recently shocked. 



3.4. Geometric Structure of Ejecta 

We observe that the ejecta plotted in the previous section appear to have a distinct 
shell-like geometry in the Southwest and Northeast regions. We attempt to characterize the 
shape of this shell with an ellipsoid by fitting bright emission with a total flux at least 15% 
that of the brightest Doppler components from all 3 regions. 

We determined the best-fit ellipsoid characterized by 8 components - the 3 spatial co- 
ordinates for the center of the ellipsoid, the 3 axes lengths, and 2 rotation angles. We 
minimized the intensity weighted RMS residuals in the 3D space by iteratively stepping 
through all plausible combinations of parameters. We show our best fit ellipsoid in Figure 
HU1 The lengths of the axes in the plane of the sky are 103"and 98.3", and the length of 
the axis perpendicular to the plane of the sky is 97.2". The ellipsoid is a sphere to within 
7%. The average residual from the best fit ellipsoid is 270 km s -1 , which is roughly 5% of 
the total velocity for ejecta on the ellipsoid. The center of this sphere is offset from the 
geometric center of the ejecta by 810 km s -1 along our line of sight. This offset can also be 
seen in the Velocity vs. Radius plot shown in Figure |3j 
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3.5. Corrugation 

Although the brightest ejecta lie near a spherical surface, the ejecta appear to be cor- 
rugated about that surface. In fact, the residuals from our best fit surface are dominated 
by systematic ~250 km s _1 corrugations about the surface rather than random small-scale 
fluctuations. This is most clearly seen in the Southwest region, as shown in Figure [TTJ We 
find that the average wavelength of the corrugation is ~24" and the amplitude is ~8" about 
the best fit surface. We further address the issue of corrugation in the next sections by 
looking at radial plot of the net intensity of the ejecta. 



3.6. Separation of Nucleosynthetic Layers 

We plot the the brightest [Si II] and [O IV] ejecta, which come from different nucle- 
osynthetic layers, for the Southwest and Northeast in Figures H] and [5j There is clearly some 
separation between these layers in some directions. We show a closeup of one such region in 
Figure [T3~| where the layers are separated by a few hundred km s _1 (corresponding to ~5"or 
~0.1 pc). The location of this region is shown in Figure ITUl 

In the SE, we observe a clump of O rich ejecta with almost no corresponding Si emission. 
We address this interesting region in more depth in §4.2. 

We further examine and quantify the separation between layers by plotting the intensity 
of the emission as a function of three-dimensional radius for all intensities of ejecta for both 
O and Si in the Southwest. Since Doppler velocity and spatial coordinates are equivalent and 
we want to determine the velocity separation between nucleosynthetic layers, we converted 
all our spatial coordinates into velocity units in order to determine the 3 dimensional velocity 
from the center of expansion fo und in §3.3. Th i s was accomplished by using the arcseconds 



to km s 1 ratio determined by iDeLaney et al.l (j2010l ). We then plotted the line flux as a 



function of 3D velocity by binning the emission in 200 km s 1 increments. 

Because we expect the behavior of the layers to vary as a function of direction, we 
created plots for many different lines of sight. For each plot we only plotted emission from a 
solid angle 7r/12 steradians wide. Our initial beam was centered 7r/8 radians above the plane 
of the sky and is wide enough to include the entire width of the region, and then incremented 
by 7r/12 radians downward for subsequent lines of sight. 

We show the flux vs. 3D velocity plots for the Southwest in Figure dD The average 
radial distance of the O with respect to Si varies significantly between locations. Along some 
lines of sight, they overlap to within one 2" bin. In different lines of sight, the O peak is at 
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a velocity up to ~500 km s 1 greater than the Si. And along yet other lines of sight, the Si 
and O peak at roughly the same velocity, but much of the O still is at larger radii than the 
Si. 

Furthermore, the peak velocity of both the Si and O chang function of direction 

from ~4400 km s _1 to ~5200 km s _1 . This is consistent with corrugation in the velocities 
of the ejecta along different directions. 

The velocity width averaged over the solid angle of the Si and O ejecta in the previous 
section are physical widths, and not just instrumental effects. We find that, to first order, 
the FWHMs are all ~1000 km s" 1 for both Si and O. The velocity uncertainty in our bins 
is roughly 130 km s _1 (since we have an uncertainty of roughly 65 km s _1 in each direction) 
and the bins are 200 km s _1 wide. 

We note that the velocity width of the peaks is much larger than the velocity width for 
any given line of sight determined in §3.1. The velocity widths of the peaks is dominated 
by averaging the 3D velocities over the entire solid angle and is not necessarily indicative of 
the velocity spread over any single line of sight. 

3.7. Faint Ejecta 

Although we dealt mostly with bright ejecta that are found to lie on a distinct spherical 
structure in the previous sections, there are weaker ejecta which lie interior to this bright 
shell, especially in the NE region. On average, these ejecta are ~10% as bright as the ejecta 
on the shock. We plot these ejecta along with the bright material in Figure [T2J We note 
that nearly all the dim ejecta lie interior to the bright ejecta. 

4. Discussion 

4.1. Supernova Explosion Physics 

The nature of core-collapse supernova explosions is a major area of research. The 
assumed structure of the star before the supernova explosion is similar for many different 
models. As a massive star fuses different elements during hydrostatic burning, it should 
produce denser and denser concentric nucleosynthesis layers, forming the classic "onion-skin" 
model of the star. 

However, between models, there is substantial variation in the relevant physics be- 
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hind the supernova explosion itself. Most groups propose neutrino-driven shocks as the 
main mechanism causing the explosion, but som e utilize diffusive, magnetic buoyancy, or 
neutrino-bubble instabilities (IJanka et al. 1 120071 ). Other groups propose jet driven explo- 
sions , where the explosion is dominated by MHD driven jets formed in rapidly rotating stars 



le.g. 



Burrows et al. 1 120071 ). 



3D maps of different nucleosynthetic layers in Cas A provide a unique opportunity to 
test and constrain the various models. Our observations focus on the layers that were initially 
near the core of the star, where the supernova explosion begins - the Fe/Ni core, the Si/S 
layer immediately above the core, and the O/Ne layer above the Si/S. We can observe the 
post-explosion geometry and velocity profile of the differ ent layers. Some mode ls predict that 
Si and O will be ejected at nearly the same velocity (e.g. iKifonidis et al. Il2006[) . while others 



predi ct that they will be ejected at velocities that differ by >500 km s 1 (e.g. IJoggerst et al. 



2009|). 



Along many lines of sight, our results are consisten t with the models of IKifonidis et al. 



( 120061 ) and the 25M models of IJoggerst et al. I ( 120091 ) - we see little difference between the 
peak velocities of O and Si. However, along other lines of sight, we find that the O and 
Si peaks are offset by hundreds of km s -1 which is inconsistent with abo ve models, but is 
consistent with the results of the 15M models of IJoggerst et al. I ( 120091 ). However, even 



though this 15M model produces a reasonable separation between O and Si, it predicts 
that the velocity width of O will be nearly twice that of Si, while we observe that both layers 
have roughly the same width along all lines of sight. Furthermore, the model predicts overall 
velocities of <2000 km s _1 , which is only half of what we observe. Put together, available 
models can reproduce most of the various behaviors that we observe, but no single model, 
so far, can reproduce observed velocity structure of the nucleosynthetic layers of the Cas A 
supernova explosion. 



4.2. The Southeast 



Ennis et al.l ( 120061 ) found a Neon "crescent" of ejecta in the Southeast and Northeast of 



Cas A where Neon was seen in the IR, but little Si was detected in the IR or X-ray. Our 
Southeast region overlaps slightly with the Neon crescent. Most of the IR emission in this 
region consists of overlapping Si and O immediately behind the X-ray ejecta. However, in 
the Neon crescent, we see only O emission that is more than an ord er of magnitude str onger 
than Si emi s sion as shown in Figure [9j This is consistent with lEnnis et al.l (120061 ) and 
Smith et al.l (120091 ) since Neon and O come from the same nucleosynthetic layer. The most 
likely explanation for this Neon/O crescent is that the Neon and O are currently encountering 
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the reverse shock and therefore becoming brighter as they are compressed, while the Si has 
yet to encounter the reverse shock. This is consistent with the picture that there is no IR 
emission from the Si because the densest Si clumps have not yet been radiatively shocked, 
and there is no observable X-ray emission because the Si ejecta have not been non-radiatively 
shocked. 



4.3. Geometry 



We observe that the center of expansion is offset from the center of the remnant by ~810 
km/s along o ur line of sight. T his is consist ent with previous results from 3D reconst ructions 
in the optical (IReed et al. Il995h and the IR ( iDeLaney et al.ll2010l ; llsensee et al. There 
are three major possible sources for the cause of this offset - asymmetries in the circumstellar 
environment, movement by the progenitor star, and the supernova explosion itself. 



Reed et al.l (119951 ) spec ulated that this offset was due to asymmetries in the pre-supernova 



circumstellar environment, llsensee et al. I (120101 ). using one patch toward the center of the 
remnant, argue that this cannot be the case since the ejecta interior to the reverse shock 
show the same velocity asymmetry, despite the fact that they are unaffected by the circum- 
stellar material. The interior ejecta are expanding into a bubble that has been cleared of 
any circumstellar material by the shocks associated with the supernov a explosion. Since ou r 
new results, which span the Bright Ring, are consistent with those of llsensee et al. I ( 120101 ). 
the velocity offset is not specific to the central regions of Cas A, but applies to the entire 
remnant. 

In principle, the velocity offset could be caused by progenitor motion, but the observed 
810 km s _1 offset is much too large a velo city for a star. Neutron stars may have velocities 



>500 km s 1 (e.g. ISatterfield et al. 1120111 ). but such large velocities are the result of "kicks" 



during the supernova explosion, w hich are caused by asymmetries in the sup ernova explosion 



itself or a binary companion (e.g. IShklovskii Ill969t IChatterjee et al. 1120051 ) 



Therefore, the observed offset it most likely caused by asymmetries in the supernova 
itself. A l ikely culprit is asymm etries formed in the first ~100 milliseconds as seen in the 



models of iBurrows et al. I ( 120071 ) and the SASI models of iBlondin et al. I ( 120031 ). Both of 
these instabilities allow the initially spherically symmetric forward shock to become highly 
asymmetric in just a few crossing times. These instabilities arise due to the response of the 
post shock pressure to changes in the shock radius. If the pressure in one region becomes 
slightly higher than the surroundings, it will push the shock outward. The preshock pressure 
drops with increasing radius, which leads to smaller pressures behind the forward shock due 
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to the outward shock displacement. If the postshock pressure radial profile is steeper than 
the preshock pressure profile, a standing acoustic wave is produced by the positive feedback 
l oop. Ejecta can "slo sh" between the standing shocks, resulting in substantial asymmetries 



( IBlondin et al. II2003I ). Presumably, the ejecta will maintain this asymmetry as they expand 



outward, resulting in a low-order asymmetry that is not necessarily centered on the location 
of the progenitor star. Note that althou gh the shock may initially be asymm etric, it will 
gradually become spherical over time (e.g. iBisnovatyi-Kogan &: Blinnikov Ill982l ). producing 
the nearly spherical surface that we observe today. 



4.4. Faint Ejecta 



The faint ejecta seen interior to the bright ejecta on the nearly spherical shell are likely 
ejecta that have yet to encounter the reverse shock. We expect to observe both [Si II] and 
[O IV] even if they are not yet shocked because they will be photoionized by energetic UV and 



X-ray photons from the reverse shock. Such ejec ta were previou s 



the remnant and in t erior to the reverse shock by lDeLaney et al.l (120101 ). ISmith et al.l ( 120091 ) 



y det e cted in the c e nter o f 



and llsensee et al. I ( I20101 ). The brightness of the ejecta appears to be roughly an order of 
magnitude less than that of nearby material which has been shocked. This is consistent with 
what is expected from a strong shock - a compression factor of about 4 is expected for a 
classic, strong, non-radiative shock, which would cause a rise in emissivity of a factor of 16. 

We note t hat the interio r eject a in the center of the remnant (see Figure [I]) discussed 
extensively in llsensee et al. I ( 120101 ) are much brighter than the interior ejecta in all the 
regions near the reverse shock. The most obvious explanation for this is that the central 
ejecta are at a higher density. This difference in density is probably caused by geometric 
effects - the central material from the remnant is traveling at about half the velocity of the 



mate rial currently encountering the reverse shock in the plane of the sky (IDeLaney et al. 



20101 ). If the ejecta are expanding homologously and were initially ejected at approximately 
the same density, the material in the place of the sky would be at a density one quarter 
that of the interior ejecta since the shocked ejecta are at twice the radius. Therefore, we 
would expect the emissivity of the central ejecta to be roughly 16 times that of material in 
the plane of the sky that is interior to the reverse shock since the emissivity of the ejecta 
varies as the square of the density. This is roughly what we observe - interior material in 
these regions have an average brightness of ~300 MJy sr" 1 , while m aterial in the center of 
the remnant have brightnesses around 4000 MJy sr _1 ( llsensee et al. 



20101 ) 
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4.5. Corrugation 



Corrugation - that is, ripples in the geometric structure of t he eject a - has been pre vious 



observed in several sup ernova remnants, including SN1006 (e.g. IWinkler fc Long 1119971 ) and 



the Cygnus Loop (e.g. iRaymond 1 120031 ). There are several possible explanations for this 



corr ugation. If the sho ck is radiative, the ripples c ould be caused b y the thermal instability 
(e.g. lBerschinger 1119861 ) or the thin shell instability ( iVishniac 1119831 ). The thermal instability 
is especially relevant for high speed (>150 km s -1 ), high temperature (T>10 5 K) shocks 
where the sound speed crossing time greatly exceeds the cooling time. This "instability" 
is actually an overstability that results from the high radiative cooling rate. The thin shell 
instability is another overstability. In this scenario, the ram pressure and thermal pressure 
are misaligned, causing ripples in the initially smooth distribution of ejecta. 

Another possibility is that the ripples are caused by inhomogeneities in the ejecta en- 
countering the shock. This inhomogeneity could be caused by Rayleigh- Taylor filaments 
created at the contact discontinuity between the ejecta and the ISM as the remnant enters 
its Sedov- Taylor phase (IWang fc Chevalier 1120011). While this is a likely explanation for the 
corrugation observed in SN1006 ( jLong et al. 1 120031 ) . it is not likely in Cas A since the ISM 
has already been swept away by the forward shock by time the ejecta encounter the reverse 
shock. 

A fin al explanation is that the corrugation is caused by density variations in the ejecta 



itself (e.g. lRaymond 1120031 ). In this model, dens e ejecta encounter ing the reverse shock would 



be slowed less than less dense clumps of ejecta. IRaymond I ( 120031 ) argues that this instability 
is the most likely cause for the corrugation in the Cygnus Loop, where the density variations 
may have been caused by ISM turbulence. 

Regardless of the mechanism, there is also a question of where the corrugation occurs. 
The observed ejecta have encountered several shocks - the forward shock during the supernova 
explosion itself, a reverse shock du ring the supernova, a nd most recently, a larger scale reverse 
shock in the supernova remnant ( jlsensee et al. 1 120101 ) . Which shock encounter creates the 
observed corrugation? 

Most studies expect corrugation at the second, larger forward or reverse shocks associ- 
ated with the remnant. However, we conclude that the most likely location for the corruga- 
tion is during one of the two shocks that the ejecta encounter during the supernova explosion 
itself. Previous studies find little evidence for decelera tion of IR ejecta after the explosion, 
including recently shocked ejecta (IDeLaney et al.ll2010l ). Our data are also inconsistent with 
recent deceleration of the ejecta. If the ejecta were suddenly decelerated, the ejecta would 
be closer to the plane of they sky in our models due to the reduction in Doppler velocity, 
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forming a flattened surface. However, we do not observe this pattern in the shocked ejecta 
- we still see a nearly perfect sphere. In other words, if the corrugation were being caused 
at the reverse shock, the ejecta would need to be decelerated by at least several hundred 
kilometers per second in order to create the observed corrugation since being shocked. This 
effect would be easily visible with our Doppler reconstruction of the global geometry, but 
is not observed. Thus, the corrugation must occur during the supernova explosion itself, 
perhaps at the initial forward or explosion reverse shock. 

5. Conclusions 

We create a 3D model of shocked ejecta of Cas A in select regions at unprecedented 
spectral resolution using IR ionic lines. We confirm previous studies that indicate that the 
remnant is offset by ~800 km s _1 along our line of sight. We find evidence for velocity 
separation between the O and Si layers along some, but not all, lines of sight. We measure 
the velocity width of these layers roughly 250 km s" 1 thick for a single line of sight, although 
the ejecta are often in bands that are ~1000 km s _1 thick averaged over several nearby lines 
of sight due to corrugation. We find evidence for corrugation in some regions of the rem- 
nant, and speculate that the corrugation was caused during the explosion itself rather than 
hundreds of years later. We use our observations of Si and O velocities to begin constraining 
models of supernova explosions, and to motivate future models to explore velocity profiles 
as a function of azimuth. 

We look forward to potential similar data sets from instruments such as the Herschel 
Space Observatory and the Stratospheric Observatory for Infrared Astronomy (SOFIA). 
Both these observatories will have the ability to create spectral cubes of Cas A at much 
higher spectral resolution. The current instruments on both observatories do not have the 
necessary instantaneous bandwidth to observe the 10,000 km s -1 velocity range of ejecta in 
Cas A, but future spectrographs will have scanning modes that will allow observations of 
ejecta at many different velocities. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. This work was supported in part by NASA/SAO Award No. AR5-6008X and 
NASA/JPL through award 1265552 to the University of Minnesota. 

K. I. would like to thank Alexander Heger for valuable insight into the physics of su- 
pernova explosions. 
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Fig. 1. — 34.81/mi [Si II] Spitzer IRS map (left) and X-ray Si Chandra map (right) of Cas A. Both 
maps have been continuum subtracted. The regions of high resolution data discussed in this text 
are indicated by the boxes. The planes shown in Figure [9] are indicated by the straight lines in the 
X-ray image. 




Fig. 2. — Typical spectra from the SH and LH Spitzer IRS module of emission in Cas A. The 
small bump near 23/im may be from the 22.9/xm [Fe III] line. 
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Fig. 3. — Velocity vs radius plot for Si (filled squares) and O (open triangle). Our assumption 
of no deceleration is good since all the ejecta lie nearly on the same semi-circle. Note that the 
semi-circle is offset from velocity by 810 km s _1 as indicated by the solid horizontal line. 
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Fig. 4. — 3D plot of the the 25.9/im [O IV] line (blue) and 34.8/zm [Si II] line (green) in the 
Southwest of the remnant as viewed from three different angles. The angles were chosen to best 
highlight the 3D structure of the ejecta. The units on the axes are arcseconds from the center of 
expansion of the remnant. 
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Fig. 5. — 3D plot of the the 25.9/mi [O IV] line (blue) and 34.8/xm [Si II] line (green) in the 
Northeast of the remnant as viewed from two different angles. The angles were chosen to best 
highlight the 3D structure of the ejecta. The units on the axes are arcseconds from the center of 
expansion of the remnant. 



Fig. 6. — 3D plot of the the 25.9/im [O IV] line (blue) and 34.8^m [Si II] line (green) in the 
Southeast of the remnant as viewed from two different angles. The angles were chosen to best 
highlight the 3D structure of the ejecta. The units on the axes are arcseconds from the center of 
expansion of the remnant. 
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Fig. 7.— 3D plot of the the 25.9/xm [O IV] line (blue), the 34.8/im [Si II] line (green), and the 
17.9/im [Fe II] line (purple) in the Southwest of the remnant as viewed from three different angles. 
The velocity axis has been stretched by a factor of approximately 1.8 in order to better highlight 
the Doppler structure of the region. The Fe emission lies on the same shell as the O and Si emission. 
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Fig. 8. — Velocity plot for the [Si II] line (dashed) over-plotted with the 25.9/xm line (solid) shifted 
under the assumption that it is either all [O IV] (left) or [Fe II] (right). The lines have been 
normalized such that the integrated flux is equal for both lines. The velocity structure matches 
very well for the assumption that the 25.9/im line is all O, but matches very poorly under the 
assumption that it is composed of Fe. 
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Fig. 9 — 3D plot of the the 25.9/im [O IV] line (blue) and 34.8/zm [Si II] line (green) in the 
Northeast (left) and Southeast (right) of the remnant along with planes of X-ray emission. The 
RA and Dec of the X-ray planes were extracted from the lines on the Chandra images in Figure [TJ 
The IR ejecta are just behind the front edge of the X-ray shock in the Northeast. The same is true 
in the Southeast, and we see a substantial amount of Oxygen dominated material where our field 
of view overlaps with the Ne-crescent detected bv lEnnis et al.l (|2006l ). This is not surprising given 
that Ne and O come from the same nucleosynthetic layer. 
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Fig. 10. — 3D plot of ejecta from all regions and the best fit ellipsoid. The units of the axes are 
km s _1 . The ellipsoid is has an eccentricity of 1.07. The zoomed in region in Figure [13]is indicated. 
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Fig. 11. — 3D plot of the 25.9//m [O IV] line in the Southwest of the remnant. The units of the 
axes are km s _1 . Two ellipsoids with radii 250 km s _1 greater than and less than the best fit 
ellipsoid are also plotted. Although the shell in this region is only ~200 km s _1 thick along any 
single line of sight, the overall velocity of the components systematically varies by ~250 km s _1 
both above and below the best fit ellipsoid. 
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Fig. 12. — 3D plot of all ejecta in the Northeast region with at least 10% the flux of the brightest 
Doppler component (left) and all ejecta with at least 25% the flux of the brightest Doppler com- 
ponent (right) . The units on the axes are arcseconds from the center of expansion of the remnant. 
The dimmer ejecta lie interior to the bright shell formed by the brightest ejecta, indicating that it 
may not yet have encountered the reverse shock. 
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Fig. 13. — 3D plot of the the 25.9fj.rn [O IV] line (blue) and 34.8/Ltm [Si II] line (green) in a select 
region of the Northeast whose location is indicated in Figure [TUl The units of the axes are km s _1 . 
The red arrow points to the center of the remnant. We detect clear separation between the O and 
Si layers along this line of sight. 
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Fig. 14. — Flux vs radius plot for the Southwest region. The O (dashed) and Si (solid) distributions 
overlap along some lines of sight, but are offset by up to ~500 km s~ x along others. The peak 
velocity of the O line and the Si line also vary from about 4400 km s _1 to 5200 km s _1 depending 
on the line of sight. This is a clear signature of corrugation. The full-widthdialf-max of the 
distributions is ~1000 km s~ x along all lines of sight. 



